THIRTY-TWO GOLDBACH VARIATIONS 



JONATHAN M. BORWEIN AND DAVID M. BRADLEY 



Abstract. We give thirty-two diverse proofs of a small mathematical gem — the funda- 
mental Euler sum identity 

C(2,1)=C(3) = 8C(2,1). 
We also discuss various generalizations for multiple harmonic (Euler) sums and some of 
their many connections, thereby illustrating both the wide variety of techniques fruitfully 
used to study such sums and the attraction of their study. 
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1. Introduction 



There are several ways to introduce and make attractive a new or unfamiliar subject. 
We choose to do so by emulating Glen Gould's passion for Bach's Goldberg variations. 
We shall illustrate most of the techniques used to study Euler sums by focusing almost 
entirely on the identities of (jl.2|) and (jl.5|) 



and some of their many generalizations. 

1.1. Euler, Goldbach and the birth of £. What follows is a transcription of corre- 
spondence between Euler and Goldbach jlH] that led to the origin of the zeta-function 
and multi-zeta values, see also [HI HE H2] • 

59. Goldbach an Euler, Moskau, 24. Dez. 1742. 1 [. . .]Als ich 
neulich die vermeinten summas der beiden letzteren serierum in meinem 
vorigen Schreiben wieder betrachtet, habe ich alsofort wahrgenommen, daft sel- 
bige aus einem blofiem Schreibfehler entstanden, von welchem es aber in 
der Tat heifiet: Si non errasset, fecerat ille minus. 2 

This is the letter in which Goldbach precisely formulates the series which sparked Euler's 
further investigations into what would become the zeta-function. These investigations 
were apparently due to a serendipitous mistake. The above translates as follows: 

When I recently considered further the indicated sums of the last two series 
in my previous letter, I realized immediately that the same series arose due 
to a mere writing error, from which indeed the saying goes, "Had one not 
erred, one would have achieved less. " 3 

Goldbach continues 

Ich halte dafiir, daft es ein problema problematum ist, die summam huius: 



in den casibus zu finden, wo m et n nicht numeri integri pares et sibi 
aequales sind, doch gibt es casus, da die summa angegeben werden kann, 



AAL: F.136, Op. 2, Nr.8, Blatt 54-55. 
2 Frei zitiert nach Marcus Valerius Martialis, I, 21,9. 
3 Opera Omnia, vol. IVA4, Birkhauser Verlag. 
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exempli grfatia], si m — 1, n — 3, denn es ist 

1 / l\ 1 / 1 1\ 1 / 1 1 1\ 7T 4 

1 + ^( 1 + 2j + ^( 1+ 2 + 3j + ^( 1 + 2 + 3 + 4j +ete - = 7^- 

1.2. The Modern Language of Euler Sums. For positive integers si, . . . , s m and signs 
<jj = ±1, consider jHj the m-fold Euler sum 

m kj 
(j . J 

((si,...,s m ;ai,...,a m ) := II 

fcl> — >fem>0 j = l ■? 

As is now customary, we combine strings of exponents and signs by replacing Sj by 7ij in 
the argument list if and only if <jj = —1, and denote n repetitions of a substring S by 
{S} n . Thus, for example, C(l) = - log 2, C({2} 3 ) = C(2,2,2) = tt 6 /7! and 

rn 

ki>--->k m >0 j=l 

The identity 

C(2, 1) = C(3) (1-2) 

goes back to Euler jHj jlHl p. 228] and has since been repeatedly rediscovered, (see, 
e.g., [2HHMlliniinZ|)- In this language Goldbach had found 

C(3,l)+C(4) = ^. 

The more general formula 

m-2 

2C(m, 1) = mC(m + 1) - ^ CO" + 1)C(™ 2 < m e Z (1.3) 

3=1 

is also due to Euler 0U ^ p. 266]. Nielsen |HI p. 229] [HI p. 198] [HI pp. 47-49] 
developed a method for obtaining (jl.3|) and related results based on partial fractions. 
Formula (jl.3|) has also been rediscovered many times [7H1 EH1 E3 HH1 1221 US] • Crandall 
and Buhler [H]] deduced (jl.3j) from their general infinite series formula which expresses 
C(s, t) for real s > 1 and i > 1 in terms of Riemann zeta values. 

Study of the the multiple zeta function (jl.lj) led to the discovery of a new generalization 
of (J1.2)) . involving nested sums of arbitrary depth: 

C({2,in = C({3} B ), neZ + . (1.4) 

Although numerous proofs of (J1.2)) and (jl.3|) are known (we give many in the sequel), the 
only proof of (|1.4|) of which we are aware involves making a simple change of variable in 
a multiple iterated integral (see [TU [THl HH| and ()5.11|) below). 
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An alternating version of ()1.2j) is 

8C(2,1) = C(3), (1.5) 

which has also resurfaced from time to time [HU p. 50] [TTJ (2.12)] [321 p. 267] and hints 
at the generalization 

8"C({2,l} n ) = C({3} n ), neZ + , (1.6) 
originally conjectured in [T3|, and which still remains open — despite abundant, even over- 
whelming, evidence jllj . 

1.3. Hilbert and Hardy Inequalities. Much of the early 20th century history — and 
philosophy — of the " 'bright' and amusing" subject of inequalities charmingly discussed 
in G.H. Hardy's retirement lecture as London Mathematical Society Secretary, [52] . He 
comments |q"2*1 p. 474] that Harald Bohr is reported to have remarked "Most analysts spend 
half their time hunting through the literature for inequalities they want to use, but cannot 
prove. " 

Central to Hardy's essay are: 

Theorem 1. (Hilbert) For non-negative sequences (a n ) and (b n ), not both zero, and for 
1 < Vi <? < 00 with 1/p + 1/q = 1 one has 

OO OO , / 

n=l m=l 

Theorem 2. (Hardy) For a non-negative sequence (a n ) and for p > 1 



7 7 — ; <ttcsc - ||a n || p ||6 n || 9 . (1.7) 

' n + m \p I 



n=l ^ ' ^ ' n=l 

We return to these inequalities in Section El 

Hardy p. 485] remarks that his "own theorem was discovered as a by-product of 
my own attempt to find a really simple and elementary proof of Hilbert 's. " He reproduces 
Elliott's proof of (|1.8p . writing 11 it can hardly be possible to find a proof more concise or 
elegant' 1 and also "I have given nine [proofs] in a lecture in Oxford, and more have been 
found since then.' 1 (See P- 488].) 

1.4. Our Motivation and Intentions. We wish to emulate Hardy and to present proofs 
that are either elementary, bright and amusing, concise or elegant — ideally all at the same 
time! In doing so we note that: 

(1) C(3), while provably irrational, is still quite mysterious, see jSHHSj and 0. Hence, 
exposing more relationships and approaches can only help. We certainly hope one 
of them will lead to a proof of conjecture (jl.6j) . 
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(2) Identities for ((3) are abundant and diverse. We give three each of which is the 
entry-point to a fascinating set: 

• Our first favourite is a binomial sum [2j that played a role in Apery's 1976 
proof, see [THIEI] and [01 Chapter 3], of the irrationality of C(3): 

C(3 ) = Ey^LJJ (1.9) 

• Our second is Broadhurst's binary B BP formula |3Uj : 

C(3) = y -7, -1, 10, -1, -7, 1, 0) + y S 3 (l, 1, -1, -2, -1, 1, 1, 0), 

where <S p (ai, a 2 , . . . , as) := Y^k=i afc2~Lp( fc+1 )/ 2 J k~ 3 , and the coefficients a fc re- 
peat modulo 8. We refer to |BJ Chapter 3] for the digit properties of such 
formulae. Explicitly, 



oo 

c(3 ) = J-Y4 



2048 11264 1024 11776 



672 2 m L(24fc+ l) 3 (24fc + 2) 3 (24A; + 3) 3 (24A; + 4) 3 
512 4096 256 3456 128 

+ 7777 777 + 7771 777 + 7777 777 + 



(24A; + 5) 3 (24A; + 6) 3 (24A; + 7) 3 (24A; + 8) 3 (24fc + 9) 3 
704 64 128 32 176 



(24& + 10) 3 (24A; + 11) 3 (24A; + 12) 3 (24A; + 13) 3 (24k + 14) 

16 216 8 64 4 

+ 



(24k + 15) 3 (24k + 16) 3 (24k + 17) 3 (24fc + 18) 3 (24k + 19) ; 
46 2 11 1 

+ 7777 7777 ~ 7777 7777 ~ 7777 7777 + 



(24& + 20) 3 (24A; + 21) 3 (24A; + 22) 3 (24fc + 23) 3 _ 

It was this discovery that lead Bailey and Crandall to their striking recent 
work on normality of BBP constants [HI Chapter 4] . 
• Our third favourite due to Ramanujan [01 p. 138] is the hyperbolic series 
approximation 

1-7 3 00 1 

C (3)- 77r 2 V 1 

^ J 180 ^ k 3 (e 2 * k - 1)' 
k=i v ' 

in which the 'error' is £(3) — 77r 3 /180 ~ —0.003742745, and which to our 
knowledge is the 'closest' one gets to writing ((3) as a rational multiple of 7r 3 . 

(3) Often results about ((3) are more precisely results about ((2, 1) or ((2, 1), as we 
shall exhibit. 

(4) Double and multiple sums are still under-studied and under-appreciated. We 
should like to partially redress that. 

(5) One can now prove these seemingly analytic facts in an entirely finitary manner 
via words over alphabets, dispensing with notions of infinity and convergence; 
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(6) Many subjects are touched upon — from computer algebra, integer relation meth- 
ods, generating functions and techniques of integration to polylogarithms, hyper- 
geometric and special functions, non-commutative rings, combinatorial algebras 
and Stirling numbers — so that most readers will find a proof worth showing in an 
undergraduate class; 

(7) For example, there has been an explosive recent interest in q-analogues, see £ 14.31 
and in quantum field theory, algebraic K-theory and knot theory, see [HJ EI] ■ 

For some of the broader issues relating to Euler sums, we refer the reader to the survey 
articles [HI EH EEl EH E3 EH EH] • Computational issues are discussed in j^J EH] and to an 
extent in [To] . 

1.5. Further Notation. 

Notation and Terminology. For positive integer N, denote the iVth partial sum of the 
harmonic series by '■= ^2n=i V n - We also use ift = T'/T to denote the logarithmic 
derivative of the Euler gamma-function (also referred to as the digamma function), and 
recall the identity ip(N + 1) + 7 = Hn, where 7 = 0.5772156649 ... is Euler's constant. 
Where convenient, we employ the Pochhammer symbol (a) n = a(a + 1) • • • (a + n — 1) for 
complex a and non-negative integer n. As usual, the Kronecker 5 mtTl is 1 if m = n and 
otherwise. 

We organize our proofs by technique, although clearly this is somewhat arbitrary as 
many proofs fit well within more than one category. Broadly their sophistication increases 
as we move through the paper. In some of the later sections the proofs become more 
schematic. We invite readers to send additional selections for our collection, a collection 
which for us has all the beauty of Blake's grain of sand 4 : 



"To see a world in a grain of sand 
And a heaven in a wild flower, 
Hold infinity in the palm of your hand 
And eternity in an hour.'''' 



2. Telescoping and Partial Fractions 



For a quick proof of (jl.2|) . consider 



n,k>0 



nk(n + k) 



1 



E^(l-^)=E^El = «3) + c(2,i) 

i,fc>0 v 7 n=l k=l 



'William Blake from Auguries of Innocence. 
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On the other hand, 

S= E \n + k) (n + k) 2 = E n (n + k) 2 + E k(n + k) 2 = 2 ^ 2, ^' 

n,fe>0 V 7 V ' n,k>0 v y n,fe>0 V ' 

by symmetry. 

The above argument goes back at least to Steinberg jH7j. See also [5%] . 
For (jl.5|) . first consider 

n=l fe=l n=l fc=l v 7 

~ti n2 ti v *(»+*) ; 



y^zn + y^z±i y 

^ nk(n + k) ^— ' n 2 (n + A;) 2 -— ' n 2 (n + fc) 

n,k>0 K ' n,k>0 v ' n,k>0 v y 



, n k/(n + k) 2 ^— ' n 2 (n + k) 

2C(2,i) + C(T,2)-C(T,2). 



Similarly, 



n=l fe=l n=l fc=l v 

oo 1 OO /■ 

n=l k=l V 



k n + k 

k fn + k- (-l) n AT 



k{n + k) 

\) k x ^ (-l) fc x ^ (-!)"+* 



y ^ + y _jz^ y 

' nk(n + k) ^ n 2 (n + k) ^ n 2 (n + k) 

n,k>0 v 7 n,fc>0 v y n,fc>0 v 7 

V + V '- 1 '"(- 1 '" +t - C (T2) 

E^P+E^p+c^-cti.s) 

n,fe>0 V ' n,k>0 v 7 

C(2,T) + C(2,T) + C(T,2)-C(T,2). 
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Adding equations (|2.1J) and (|2.2J) now gives 

2C(2,1) = C(3) + C(3), 

i.e. 

oo 



771=1 V ' 



n=l 



which is (|1.5j) . 



3. Finite Series Transformations 
For any positive integer N, we have 

AT „ AT , n-1 AT , n 



El ^ — v 1 ^ — v 1 ^ — v 1 ^ — V 1 
n 3 Z^ n 2Z^^ ^U 2 ^JV-fc+l 



n=l n=l k=l n=l k=l 

by induction. Alternatively, consider 

N „ JV 



r:= y 1 = y(l-l) 7 JL- 

' nk(k — n) ^— ' kl{k — nr 

n,k=l v 7 n,fc=l V 7 V y 



On the other hand, 

^— ' ri 2 \ k — n k 

n,h=l v 

" i i " i All 

n=l x fe=l fc=n+l k=l 

N N n-1 N ^ / N t AT . 

A^ n 3 ^ n 2 ^ n _j i ^/ n 2\ A^/ k-n ^ I 

n=l n=l k=l n=l N fc=n+l fe=l 

Since T = 0, this implies that 

AT j ^ 1 n_1 1 ^ 1 / ^ 1 Af_n 1 \ ^ 1 " 

/ r)3 < r)2 < i* < r)2 ( < t- < 1* ) < r)2 < 



n 3 2 -— ' n 2 ' ' n 2 V ' A; ' 7 ' n 2 ' N 

n=l n=l fc=l n=l N fc=l fe=l ' n=l fc=l 
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which is (|3.1jl . But the right hand side satisfies 

u N 1 N n 1 

N ~^n 2 'N-^n 2 ^N-k + l 

n=l n=l fc=l 

El n 1 ^ / I 1 

' jV-n + 1 ~ iVTI^ In + N -n + l) 

n=l n=l x ' 

Letting A" grow without bound now gives p.2jl . since lim — — = 0. 

4. Geometric Series 



N + l 



□ 



4.1. Convolution of Geometric Series. The following argument is suggested in [70] . 
A closely related derivation, in which our explicit consideration of the error term is sup- 
pressed by taking A" infinite, appears in [22] • Let 2 < m G Z, and consider 

m-2 N N m-2 

V C(j + l)C(w - j) = lim V V V — -- r 

j=l n=l fc=l j=l 

= lim I V f 1 1 \ m-2 

k^n 

N 



n m+l 



(m - 2)C(m + 1) + 2 lim ^ ) 

n,k=l v y 



Thus, we find that 



m-2 

(m - 2)C(m + 1) - J] CO' + l)C(m - j) 

3=1 

N „ JV 



, " 1 " A 1 

2 lim V — V - - 

at^oo ^ n m ^ V k k- 

n=l k=l v 

k=£n 

" 1 r^i i " i 
2 lim > — O ---+> 

n=l k fc=l fc=l 

N n 

2C(m, 1) - 2((m + 1) + 2 lim V — V 



n=l k=l 
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and hence 



m-2 



N n 



1 1 

2C(m, 1) = mC(m + 1) - £ C (j + l)C(m - j) - 2 jim £ — £ N _ k + r 

j=l n=l k=l 

But, in light of 

y^y i <y-- - <—y ( - + -) 

n=l ' 



IE 



N 



n=l k=l 



n=l 



N+V 



the identity (J1.3|) now follows. 



□ 



4.2. A Sum Formula. Equation (jl.2j) is the case n = 3 of the following result. See [25] . 
Theorem 3. J/3 < n G Z t/ien 



n-2 



(4.1) 



We discuss a generalization (|l().4jl of the sum formula (|4.1jl to arbitrary depth in §10.21 
Proof. Summing the geometric series on the right hand side gives 

n—2 oo oo _ oo 

yvy - = y 

^ h? (h + m) n -i 

=1 h=l m=l K ' 



h n 2 m(h + m) m(h + m) n 1 
h,m=l v ' v ' 

V — V (- —) -C(n-1 1) 

h=l m=l v 7 

h 



E^E 5 -««-m) 

/i=l fc=l 
x. ^ v ^ n— 1 

E^+E^E^"- 1 . 1 ) 

h=l h=l k=l 

C(n) 



□ 



4.3. A g-Analogue. The following argument is based on an idea of Zudilin (7H|. We 
begin with the finite geometric series identity 



uv 



uv 



uv 



s-1 



uv 



(l-u)(l-uv) s (l-v)(l-uv) s (l-u)(l-v) s (1 - v)i +1 (l - uv) s ~i 
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valid for all positive integers s and real u, v with u ^ 1, uv ^ 1. We now assume s > 1, 
q is real and < q < 1. Put it = g m , v = q n and sum over all positive integers m and n. 
Thus, 

„m+n „m+2n 

y ? + y q 



(1 - q m )(l - q m+n ) s ^ (1 - q n )(l - q m+n ) s 

m,n>0 K * /V * / m,n>0 v * 7V * 7 
„m+n „m+2n 

Q \ ^ q 

(1 - q m )(l - q n ) s ^ (1 - q n ) s (l - q m + n ) 

m,n>0 K * 7V H 1 m,n>0 K H > K * 7 
s ~ 2 ,m+2n 



EE 



(1 - q n V +l (l - q m + n y~i 

3=1 m,n>0 v 1 J \ 1 ) 



E 9 Q Q _ \ ^ 

m,n>0 ^ " 



s - m+2n 



q 



s-2 



j=l m,n>0 

m+2n 



fx _ guv 2-^ x — o m ^— ^-^ (l — o n V' +1 Cl — (7 m+n )' s ^' 

n>0 v ^ 7 m=l ^ 3=1 m,n>0 v 1 / V 1 / 

q2n qU+m \~~^ qm+2n 



V — + V V V — 

Z-^ (i _ Q n\s+i (1 — n n ) s (l — a m ) ^— ' (I — a n V +1 (l — a m + n ) s -i 

n>0 v ^ 7 ra>m>0 v * ' v ^ 7 j=l m,n>0 v 1 / V 1 7 



Cancelling the second double sum on the left with the corresponding double sum on the 
right and replacing m + n by k in the remaining sums now yields 

9 \^ q \ ^ \ ^ q 

(1 — q h ) s (l — q m ) ~ h — qny+i (i _ q m)j+i(i _ q ky-j ' 

fc>m>0 V y 7 V ^ 7 n>0 V H 7 j=l fc>m>0 V H ' K H ' 

or equivalent ly, that 



V - = V - + T V - (4 2) 

^(l-q k Y+ 1 ^ (l-q k ) s (l-q m ) ^ ^ (1 - g*) s -i(l - q m )3+ 1 ' K ' J 

k>0 V ^ 7 fc>m>0 V ^ 7 V ^ 7 J=l fc>m>0 V M / V * / 

Multiplying (|4.2)l through by (1 — q) s+1 and letting q — > 1 gives 

s-2 

C( s + i) = C(^i) + E^-^^' + 1 )' 

which is just a restatement of (|4.1j) . Taking s = 2 gives ()1.2|) again. □ 
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As in |21], define the g-analog of a non-negative integer n by 



n-l 



and the multiple g-zeta function 



l-q" 



l-o 

k=0 * 



c[ Su ..., Sm }:= nv?^' (4 - 3) 

where Si,S2,...,s m are real numbers with si > 1 and Sj > 1 for 2 < j < m. Then 
multiplying ()4.2|) by (1 — q) s+1 and then setting s = 2 gives ([2, 1] = £[3], which is a 
g-analog of (jl.2j) . That is, the latter may be obtained from the former by letting q — > 1— . 
On the other hand, s = 3 in ()4.2|) gives 

C[4] + (1 - q)C[3] = C[3, 1] + (1 - q)C[2, 1] + C[2, 2], 

which, in light of C[2, 1] = C[3] implies C[3, 1] = C[4] - C[2,2]. By Theorem 1 of |23|, 
we know that C[2,2] reduces to depth 1 multiple g-zeta values. Indeed, by the g-stume 
multiplication rule [23, C[2]C[2] = 2([2,2] + £[4] + (1 - g)([3]. Thus, 

C[3, 1] = C[4] - C[2, 2] = |C[4] - \ (C[2]) 2 + |(1 - g)C[3], 
which is a g-analog of the evaluation [HI H3 CHI HH1 C3 1201 12H 

((3,1) = — . 
sv ; 360 

Additional material concerning g-analogs of multiple harmonic sums and multiple zeta 
values can be found in W\ I2B | l27j . 



5. Integral Representations 



5.1. Single Integrals I. We use the fact that 

= k 2 



1 1 



u k - 1 (-\ogu)du = k>0. (5.1) 



Thus 

fc>n>l n=l k>n u 

= E- / i-^gu)^- 1 du 

n=l J ° k>n 



f^/ ( - l0gn) l^ M 
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- / > — du 

/ n 1 — it ^— ' n 

J0 n=l 

/ (— logu)(l — u)^ 1 log(l — w) _1 dn. 
Jo 



(5.2) 



The interchanges of summation and integration are in each case justified by Lebesgue's 
monotone convergence theorem. After making the change of variable t = 1 — u, we obtain 



S >n>l J " /0 n=l 



^n— 1 

d*. (5.3) 

Again, since all terms of the series are positive, Lebesgue's monotone convergence theorem 
permits us to interchange the order of summation and integration. Thus, invoking (|5.1j) 
again, we obtain 

£ ^ = E i / (-logo^^E-T' 

fc>n>l n=l ^ u n=l 

which is (JH2I). □ 
5.2. Single Integrals II. The Laplace transform 

1 '°° r((j + i) 



x r - l {-\ogxf dx= e~ ru u a du = K ! , r > 0, a > -1, (5.4) 
n Jo r ' a 



generalizes ()5.1|) and yields the representation 

T(m+1) _ 1 ^ ri , . mj _ (-1)- ^log" 1 * 



i ^r(m+i) i ^ r i r _ lf (-ir r 

C(m+1) = — - > — — = — - > / x r L (-\ogx) m dx = - — 

ml z — ' r m+i m z — / I m / 

r=l r=l 1/0 J " 



1 — X 



The interchange of summation and integration is valid if m > 0. The change of variable 
x i — > 1 — x now yields 

£( m + l) = ± > / log^M-a:)—, l<meZ. (5.5) 

m! y x 

In ^H], equation (|5.4|) in conjunction with clever use of change of variable and integra- 
tion by parts, is used to prove the identity 

oo oo oo 1 

l! « lt! )=EI-E^ £ ? . »^^. (5.6) 

"1=1 "2=1 "fc=l p=l+rti+n 2 H h«fe 

The case = 1 of (|5.6j) is precisely (jl.2|) . We give here a slightly simpler proof of (|5.6|) . 
dispensing with the integration by parts. 
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From 1)5. 4|L 

«C(* + 2) = gi.^ = gi J [V'(-lcg 

= / (- logx) fc log(l - x) 
Jo 



x) fc cfe 



_ x <ix 

a; 

1 dx 

W^l - X^f-logx) 



1 X 

oo oo oo „ ,! ^ ni+n2+ ... +nfc 



rai=l n2=l n 
oo oo oo 



1 


nin 2 ■ ■ 


■ ^fe 


1 




nin 2 ■ ■ 


•n* 


1 




nin 2 • ■ 





/ x p ~ 1 {—\ogx)dx 

-n 2 +-+n k ^° 

E ^ 



ni=l 712=1 "fc=l p>ni+n2H hn fc 

oo oo oo 



"1=1 U2=l nfc=l p>ni+n 2 H \-n k J 



□ 



5.3. Double Integrals I. Write 

E s^-n^^EM™- 1 ^ 

fc,m>0 v y fc>0 m>0 

Jo Jo l-xy 

Now make the change of variable u = xy, v = x/y with Jacobian l/(2v), obtaining 

2 J 1 - u J u v J l-u 

which is (|5.2|) . Now continue as in ^5.11 □ 

5.4. Double Integrals II. The following is reconstructed from a phone conversation 
with Krishna Alladi. See also j^j. Let e > 0. By expanding the integrand as a geometric 
series, one sees that 

1 f f (xy)' 



n=l 

Differentiating with respect to e and then letting e = gives 

" J \og(xy) ^ 1 /" /" logx + logy j j f f 1 \ogx ^ 



c(3) = -i /' r^) iI% = -i r r^iy dxdy = -rr 

2 J J l-xy 2 Jo Jo l-xy Jo Jo 



1 — xy 
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by symmetry. Now integrate with respect to y to get 

C 1 dr 
C(3)= / (logx)log(l-x)-. (5.7) 
Jo x 

Comparing (|5.7|) with ()5.3|) completes the proof of ()1.2|) . □ 

5.5. Integration by Parts. Start with ()5.7|) and integrate by parts, obtaining 

2C 3 = / -^dx = / log 2 (l-x)— = > / — dx= > — -. 

w Jn 1-x Jo x ^ Jn nk ^ nk(n + k) 

J0 JU n,k>0 Jl> n,k>0 V ' 

Now see §2 □ 

5.6. Triple Integrals I. This time, instead of (j5.1|) we use the elementary identity 

1 _ / „,-i / 
k 2 n 

This yields 



/ yf 1 / v\~ n ~ X \ Vl' 1 dy 3 dy 2 d yi , k>n>0. 
Jo Jo Jo 



1 /"l /"l/i /"2/2 

Y,VT= V? (i-y^r 1 {l-y,)- l dy,dy 2 d yi . 

r> K 71 JO JO JO 



(5i 



fc>n>0 



Now make the change of variable yi = 1 — Xj for i = 1, 2, 3 to obtain 

To - = / (1 — / Xn 1 [ x^ 1 dx 3 dx 2 dxi 

k 2 n Jr. J Tt i T „ 

= / Xg 1 / x 2 x / (1 — dx\ dx 2 dxz- 
Jo Jo Jo 

After expanding (1 — xi) -1 into a geometric series and interchanging the order of sum- 
mation and integration, one arrives at 

k^ = J2 x 2 1 / ^r 1 ^2 = X) --3 ' 

fe>n>0 n=l J ° J ° J ° n=l 

as required. □ 
More generally [HI Q31 [TBI EEH1 HB ES] , 



c(-x,...,- fc )= e IK /n(n 

ni>--->n fe >0 j=l j=l ^ r=l 



(5.9) 



where the integral is over the simplex 

i>t«>.-.>4;)>.--># ) >--->C>o, 
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and is abbreviated by 

The change of variable 1 1— > 1 — £ at each level of integration switches the differential forms 
a and 6, thus yielding the duality formula [Hj [23 P- 483] (conjectured in [33]) 

C(si + 2, . . . , Sn + 2, = C(r n + 2, {1}% . . . , n + 2, {l}* 1 ), (5.11) 

which is valid for all nonnegative integers si, ri, . . . , s n , r n . The case si = 0, r\ = 1 
of (|5.11|) is (jl.2j) . More generally, ()1.4|) can be restated as 

I (ab 2 ) n = [ {a 2 b) n 
Jo Jo 

and thus (|l-4j) is recovered by taking each Sj = and each t$ = 1 in (|5.11|) . For further 
generalizations and extensions of duality, see [TSJ |2H ■ 

For alternations, we require in addition the differential form c := —dt/(l + t) with which 
we may form the generating function 

^^ 3 "C({2,1}") = Y J \z 6n+3 / (ac 2 ab 2 ) n ac 2 + z 6n+6 / (ac 2 ab 2 ) 6n+e \. 

n=l n=0 I JO Jo J 

A lengthy calculation verifies that the only changes of variable that preserve the unit 
interval and send the non-commutative polynomial ring Q(a, b) into Q(a, b, c) are 





b) 


= S(a,b), 


i i-> 




(5.12) 


S(a. 


b) 


= R(b, a), 




1-t, 


(5.13) 


S(a. 


b) 


= S(2a,b + c), 


t i-> 


t 2 , 


(5.14) 


S(a. 


b) 


= S(a + c,b — c), 


t ^ 


2t 
1 + t' 


(5.15) 


S(a. 


b) 


= S(a + 2c, 2b -2c), 


t i-> 


4t 

(1 + t) 2 ' 


(5.16) 



and compositions thereof, such as t i— > 1— 2t/(l + 1) = (1— £)/(l+i), etc. In ()5.12|) - (j5.16|) . 
S(a, b) denotes a non-commutative word on the alphabet {a, b} and R(b, a) denotes the 
word formed by switching a and b and then reversing the order of the letters. 

Now view a, b and c as indeterminates. In light of the polynomial identity 

ab 2 - Sac 2 = 2[ab 2 - 2a(b + c) 2 } + 8[ab 2 - (a + c)(6 - c) 2 } + [(a + 2c) (26 - 2c) 2 - ab 2 ] 

in the non-commutative ring Z(a,b, c) and the transformations ([5.14)1 . (|5.15j) and (|5.16|) 

above, each bracketed term vanishes when we make the identifications a = dt/t, b = 



18 



J. M. BORWEIN AND D. M. BRADLEY 



dt/(l —t),c = —dt/(l + t) and perform the requisite iterated integrations. Thus, 

C(2,1)-8C(2,1)= f ab 2 -8 [ ac 2 = 0, 
Jo Jo 

which in light of (jl.2|) proves (jl.5j) . □ 

5.7. Triple Integrals II. First, note that by expanding the integrands in geometric series 
and integrating term by term, 

£/ 2 !) = 8 i* 1 dx f x y d y f v zdz 



x Jo 1-2/ io 1 - z 
Now make the change of variable 

xdx du y dy dv z dz dw 



1 — x 2 1 + u 1 — y 2 1 + v' 1 — z 2 1 + w 

to obtain the equivalent integral 

f°° ( du du du \ f u dv f v dw 

C(2, 1) = 8 / — + 



/ \2u 2(2 + u) 1 + uJJq 1+v J l + w' 
The two inner integrals can be directly performed, leading to 



J u{u+l)(u + 2) 



Finally, make the substitution u + 1 = 1/ yl — x to obtain 



c( 2,i) = i/ 1 !2i!(l^) dI = c( 3), 



by Q . 



5.8. Complex Line Integrals I. Here we apply the Mellin inversion formula [3J p. 243], 
[721 pp. 130-132 ] 

i r +iQo d Z I 1 ' y> 1 

which is valid for fixed c > 0. It follows that if c > and s — 1>c>1— t then the 
Perron-type formula 

' n\ 2 dz 



1 °° °° -r rc+ioo / 

as,t) + -as + t) = j2n-^k- t ^- / 

Z n=l fc=l Z7F< •'c-ioo V 



Z 



"C+IOD 1 

2^/ C(«"*)C(* + *)^ (5-17) 



GOLDBACH VARIATIONS 



19 



is valid. (Interchanging the order of summation and integration is permissible by absolute 
convergence.) Although we have not yet found a way to exploit (|5.17jl in proving identities 
such as fll.2j) . we note that by integrating around the rectangular contour with corners 
(±c±iM) and then letting M -> +00, one can readily establish the stuffle [T51 ITTH 1231 
formula in the form 

C(S, t) + ^((S + t) + C(t, S) + ^((t + S ) = ({s)((t), S,t>l + C. 

The right hand side arises as the residue contribution of the integrand at z — 0. One can 
also use ()5.17j) to establish 



00 _ 00 

s=2 

but this is easy to prove directly. 



x 



mn(n — x) 2 ^— ' n(n — x) 

n>m>0 y ' n=l v ' 



5.9. Complex Line Integrals II. We let X(s) := 2~2 n >o X n n~ s represent a formal Dirich- 
let series, with real coefficients A„ and we set s := a + i r with a = 3?(s) > 0, and consider 
the following integral: 

2 



POO 


X(s) 


2 -1 poo 

dr=- 


X(s) 


Jo 


s 


2 7-OC 


s 



dr. 



(5.1* 



as a function of A. We begin with a useful variant of the Mellin inversion formula, namely 

(5.19) 



COs{at) -dt=*e- 



t 2 + u 2 



for u,a > 0, as follows by contour integration, from a computer algebra system, or 
otherwise. This leads to 

Theorem 4. (Theorem 1 of [Zj). For X(s) = Y2™=i An^ _s and s = a + ir with fixed 
a = 3?(s) > such that the Dirichlet series is absolutely convergent it is true that 



POO 


X(s) 


Jo 


s 



dr 



2a 



E 

n=l 



2rr 



(5.20) 



where A n := Ylk=i and A := 0. 

More generally, for given absolutely convergent Dirichlet series a(s) := Y2n°=i a n n 
and p(s) := J2n=i ^n' 8 

a(s) f3{s 



a 2 + r 2 



dr 



71 

2a~ 



E 

71=1 



A n B n — A n _i -B n -1 



n 



2cr 



(5.21) 



in which A n = Y2=i a k and B n = Y2=i Pk- 
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Note that the righthand side of (|5.20J) is always a generalized Euler sum. □ 
For the Riemann zeta function, and for a > 1, Theorem 0] applies and yields 

- k (a) = C(2(7-l)-ic(2a), 

7T 2 

as A n = 1 and A n = n — 1/2. By contrast it is known that on the critical line 

1/2 /1\ , , / — s 1 
^ t< (-)=log(v^)--o, 

There are similar formulae for s ^ ((s — k) with k integral. For instance, applying the 
result in (fFT2T)j) with £i := t i-> ((t + 1) yields 

on using (|1.2j) . For the alternating zeta function, a := s i— > (1 — 2 1 ~ s )£(s), the same 
approach via (j5.21j) produces 



7T 

and 



1 p(3/2 + fr)«(3/2 + fr ) <fr = + = 3 _ 9 

7T Jo 1/4 + r z 4 



_L £ n (3/2 + frK(3/2 + fr) dT = 1} + C(2J) + q , (3) = 9 c(2) , og(2) _ 3 ((3)> 

since as we have seen repeatedly C(2, 1) = C(3)/8; while C(2,T) = C(3) - 3/2 C(2) log(2) 
and C(2,T) = 3/2C(2)log(2) - 13/8 C(3), (e.g., |TZj). 

As in the previous subsection we have not been able to directly obtain (jl.5j) or even (jl.2|) . 
but we have connected them to quite difficult line integrals. 

5.10. Contour Integrals and Residues. Following let C n (n G Z + ) be the square 
contour with vertices (±1 ± i)(n + 1/2). Using the asymptotic expansion 

i){z) ~logz-^- j^JpF* |arg^|<7r 

" " r=l 



in terms of the Bernoulli numbers 



1 - e~* 2 ^ (2r 

r=l v 

and the identity 



1p(z) = lf){—z) COt HZ, 
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we can show that for each integer k > 2, 



n— >oo 



lim / z~ k i)\-z)dz = §. 



Then by the residue theorem, we obtain 

Theorem 5 (Theorem 3 of jHS])- For every integer k > 2, 

oo k— 1 

2 il>(n) = k((k + 1) - 2j((k) - C0')C(* - i + 1) 

n=l jf=l 

where 7 = 0.577215664. . . is Euler's constant. 



In light of the identity 



fe=i 

Theorem is equivalent to (jl.3|) . The case k = 2 thus gives (jl.2|) . □ 

Flajolet and Salvy |UJ developed the residue approach more systematically, and applied 
it to a number of other Euler sum identities in addition to (|1.3j) . 

6. WlTTEN ZETA-FUNCTIONS 
We recall that for r,s > 1/2: 

00 00 

W(r,M) : =EE , w 7 v 
z — ' z — ' n r m s (n + mr 

n=l m=l v ' 

is a Witten ^-function, [7711121110]. We refer to |77j for a description of the uses of more 
general Witten ^-functions. Ours are also called Tornheim double sums, [Hi]. There is a 
simple algebraic relation 

W(r,s,t) =W(r-l,s,t + l) + W(r,s-l,t+l). (6.1) 

This is based on writing 

m + n m n 



(m + n) t+1 (m + n) t+1 (m + n) t+1 
Also 

W(r,s,t) = W(s,r,t), (6.2) 

and 

W(r, s, 0) = C (r) C(s) while W(r, 0, t) = r). (6.3) 
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Hence, W(s, s,t) = 2 W{s, s - 1, t + 1) and so 

W(l, 1, 1) = 2 W(l, 0, 2) = 2 C(2, 1) = 2 C(3). 

Note the analogue to l|P|) . viz. £(s, *) + C(*, s) = C(s) C(*) - C(« + *)> shows °> s ) = 
2 C(a, s) = ( 2 (s) - ((2s). Thus, W(2, 0, 2) = 2 C(2, 2) = tt 4 /36 - vr 4 /90 = vr 4 /72. 

More generally, recursive use of 1)6.1)1 and (j6.2j) . along with initial conditions (jOj) shows 
that all integer W(s, r, i) values are expressible in terms of double (and single) Euler sums. 
If we start with Y(s)/{m + n) 1 = J (— logo - )' -1 a m+n ~ 1 da we obtain 

W(r, 8,t) = ±r f'uAa) Li» ^ loga ^ _ da. (6.4) 
For example, we recover an analytic proof of 

2 C(2, 1) = W(l, 1,1)= / ^ (1 - a) da = 2 C(3), (6.5) 

Jo a 

Indeed S in the proof of §21 is precisely W(l, 1, 1). 

We may now discover analytic as opposed to algebraic relations. Integration by parts 
yields 

W(r, 8 + 1, 1) + W(r + 1, s, 1) = Li r+1 (l) Li s+1 (l) = ((r + 1) ((s + 1), (6.6) 

So, in particular, W(s + 1, s, 1) = ( 2 (s + l)/2. 

Symbolically, Maple immediately evaluates W(2, 1, 1) = 7r 4 /72, and while it fails di- 
rectly with W(l,l,2), we know it must be a multiple of 7r 4 or equivalently C(4); and 
numerically obtain W(l, 1, 2)/C(4) = .49999999999999999998 . . .. 

6.1. The Hilbert Matrix. Letting a n := l/n r and 6 n := l/n s , inequality (|1.7jl of Section 
O yields 



W(r, s, 1) < 7T esc (J J VCM. (6.7) 

Indeed, the constant in (|1.7jl is best possible [2211101 • We consider 

W((p-i)*,M) 

and observe that with a^(s) := Ylm=A n / m )~^'~ V>S I \ n + m ) ~~ * ^csc > we have 

°° °° n _s m~ (p ~ 1)s °° 1 
£ p : = lim (ps - 1) > > = lim Ops - 1) > ^(s) 

n=l m=l n=l 



oc 



. \ — VTCSC (7r/o))} ,. ,„ 1NW . / 7T 



n=l 
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+ 7T CSC 



Setting r := (jp — l)s,s — ► we check that C,{ps) l / p (^(qr) 1 ^ = C(p s ) an d hence the 

best constant in (|fi.7|) is the one given. 

To recapitulate in terms of the celebrated infinite Hilbert matrix, Ho : = {l/(m + n)}^ n=l , 
pp. 250-252], we have actually proven: 

Theorem 6. Let 1 < p, q < oo be given with l/p + 1/q — 1. T7ie Hilbert matrice Tio 
determines a bounded linear mappings from the sequence space i v to itself such that 

/to p,p = llm 77 v = 7TCSC - 

»-i/p C(ps) VP 

Proof. Appealing to the isometry between (£ p )* and £ q , and given the evaluation £ p 
above, we directly compute the operator norm of 7io as 

fir 

\\7~(-o\\p,p — sup HWo^Hp = sup sup (Hox,y) = n esc I — 

|^|| P =1 lls/||<j=l ll a; llp= 1 \P 

□ 

A delightful operator-theoretic introduction to the Hilbert matrix 7io is given by Choi 
in his Chauvenet prize winning article |37j . 

One may also study the corresponding behaviour of Hardy's inequality (jl.8j) . For 
example, setting a n := 1/n in ()1.8|) and denoting H n := Yuk=i yields 

Application of the integral test and the evaluation 

iogx\ p , r (i + p) 

' ax — 



x 



(p-1 



P+i ' 



for p > 1 easily shows the constant is again best possible. 

7. A Stirling Number Generating Function 

Following [35], we begin with the integral representation ()5.5j) of £ 15.21 In light of the 
expansion 

± ^log m (l - x) = Vu(n,m) — , 0<meZ, 

ml nl 

n=0 
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in terms of the unsigned Stirling numbers of the first kind (also referred to as the Stirling 
cycle numbers in [ED]), we have 

f 1 [ / x n> \ dx u ( n m \ 

C(m+1)= \J2 u ^ m ^ - = E A ^' l<meZ. 

Jo ^ n=l ' ' n=l 



n\ n 



Telescoping the known recurrence 

u(n, m) = u(n — 1, m — 1) + (n — l)w(n — 1, m), 1 < m < n, (7-1) 

yields 

u(n, m) = (n - 1)! j<W + E " '™"^ ! • ( 7 - 2 ) 
Iterating this gives the representation 



7 

3=1 J 



C(m+l)=C(2,{l} m - 1 ), l<mGZ, 

the m = 2 case of which is (|1.2)1 . See also n = in (jlO.HJ) below. □ 

For the alternating case, we begin by writing the recurrence (17.1)1 in the form 

u(n + 1, k) + (j — n)u(n, k) = u(n, k — 1) + j u(n, k). 

Following [S3], multiply both sides by (—\y i + k + 1 j k - rn - 1 j(^j _ ra ) n; where 1 < n < j — 1 
and k,m G Z + , yielding 

fc f (-1)" +1 M (n + l,fc) _ (-l)"«(n,fc) 1 ;fc _ m _ 1 
I (j - n) n (j - n + l) n _i j 

C-F) n 

- 1 ; {(-l) k - l u{n,k-l)j k - m - 1 - {-if u{n,k)f- m } . 



(j - n) n 

Now sum on 1 < k < m and 1 < n < j — 1, obtaining 

1 (-1)^+1 ^ 
-T-^ = (—1)1 Z^" 1 ) 0-n-l)!w(n,m). 

Finally, sum on j e Z + to obtain 

k=l j=k n=m j=n+l ^ J ' 

Noting that 

y (j-l-n)! f> k\ i pni n 1 

i£ti + n! 2U '' n+ ' J (n-l)!(n-l) : 
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we find that 



c(m) = jrjt { ~ l)3+ku{3 ' k) + jr (_1)n+m M(n ' m 



j\ jm-k L^i < n _ i)i ( n - 1) 

fe=l j=fc J J n=m K 1 v y 

Now employ the recurrence (|7.1|) again to get 
c(m) _\- 2 y- (- 1 ) J+fc ^ fc ) , | ^ (-!)'+" M (j,m) 

' o'l jm—k / j j\ j / j j\ 

k=l j=k J ' J j=m-l J ' J j=m J ' 

°^ l,m) ^ (-l) n+m u(n- l,m- 1) 

+ ^ (n-1)! + ^ (n-l)!(n-l) 

ra=m v ' n=m v y v ' 

_^^ (-iy +k u(j,k) 2 ~ (-l)^-^( J; m-l) ^ 

o'l jm—k / j j) j V • / 

fc=l j=fc J ' J j=m-l J ' J 

Using (|7.2jl again, we find that the case m = 3 of ()7.3j) gives 
„ x v 2 - u(j, 1) ^ 2) 

c(3) = E - ? , ?2 J + 2 E ? , ? 

-f:^ +a x:^-i).£sM 

i=i J j=i J ' J k=i 

_f^ (-l) i+1 , 2 f>(-l) J ^l 

i=i J i=2 J fc=i 



= 2C(2,1)-C(3), 

which easily rearranges to give (|2.3|) . shown in ^21 to be trivially equivalent to (jl.5|) . □ 

8. POLYLOGARITHM IDENTITIES 

8.1. Dilogarithm and Trilogarithm. Consider the power series 

n 2 k ' 



J(x) :=Cx(2,l)= 4r> 0<x<1. 



n>fc>0 

In light of ()12.3|) . we have 



o t J 1 - u J 1 - v J 2t 



The computer algebra package Maple readily evaluates 

f log (1 - t) dt = ^ + 1 _ ^ + _ x ^ Li2 ^ _ ^ _ Lia ^ _ ^ ^ ^ 



'0 



2f 
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where 

oo 

Li s (x) := 



n" 

n=l 



is the classical polylogarithm [STfl I60j . (One can also readily verify the identity (j8.1|) by 
differentiating both sides by hand, and then checking ()8.1J) trivially holds as x — > 0+. See 
also p. 251, Entry 9].) Thus, 

J(x) = C(3) + - log 2 (l - x) log(z) + log(l - x)Li 2 (l - x) - Li 3 (l - x). 

Letting x — > 1— gives (jl.2|) again. □ 
In p. 251, Entry 9], we also find that 

J{-z) + J(-l/z) = -| log 3 z - Lia(-z) log z + U 3 (-z) + C(3) (8.2) 

and 

J{l-z) = ilog 2 zlog(^-l)-|log 3 z-Li 2 (l/z)logz-Li 3 (l/^) + C(3). (8.3) 
Putting z = 1 in ()8.2j) and employing the well-known dilogarithm evaluation |H9*1 p. 4] 

^ -ir 7T 2 



Lia(-l) = £ 



n 2 12 



n=l 

gives (jl.5|) . Putting z = 2 in (j8.3|) and employing the dilogarithm evaluation [SHI p. 6] 

L 2 - = > — — = bg 2 2 

2 \2/ ^n 2 2 n 12 2 s 

v 7 n=l 

and the trilogarithm evaluation [HHl p. 155] 

Li 3 - = > — — = -C 3 log 2 + - bg 3 2 

3 \2 / ^ n 3 2 n 8 SV ; 12 s 6 S 

x ' n=l 

gives f)l .5j) again. □ 
Finally, as in fHJ Lemma 10.1], differentiation shows that 

j <-*> = - jw + + ivh) - Hj^w) ■ < 8 - 4 » 

Putting [T51 Theorem 10.3] x = 1 gives 8J(— 1) = J(l) immediately, i.e. (jl.5|) . □ 

In [T^j, it is noted that once the component functions in (J8.4|) are known, the coefficients 
can be deduced by computing each term to high precision with a common transcendental 
value of x and then employing a linear relations finding algorithm. We note here a 
somewhat more satisfactory method for arriving at ()8.4|) . 

First, as in £ 15 .61 one must determine the fundamental transformations (|5.12j) - (j5.16|) . 
While this is not especially difficult, as the calculations are somewhat lengthy, we do not 
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include them here. By performing these transformations on the function J(x), one finds 
that 



It now stands to reason that we should seek rational numbers n, ri, r% and r 4 such that 

ac 2 = r x ab 2 + 2r 2 a(b + c) 2 + r 3 (a + c)(b- cf + r 4 (a + 2c)4(6 - cf 

is an identity in the non-commutative polynomial ring Q(a, b, c). The problem of finding 
such rational numbers reduces to solving a finite set of linear equations. For example, 
comparing coefficients of the monomial ab 2 tells us that ri + 2r2+r^+Ar^ = 0. Coefficients 
of other monomials give us additional equations, and we readily find that r\ = — 1, 
r 2 = 1/4, r 3 = 1 and r 4 = —1/8, thus proving ()8.4j) as expected. 

8.2. Convolution of Polylogarithms. Motivated by [USES], for real < x < 1 and 
integers s and t, consider 




oo 



.n+m 



oo 



: n+m {m — n + n) 
n s m t+1 (m — n) 




x 



E 



x 



t (m — n) 



m,n=l 





U s (x)Li t+1 (x) - Li s+t+1 (x 2 ) +T s _ lji+ i(x). 



Telescoping this gives 



s 



T S)t (x) = T 0jS+ i(x) - sLi s+t+1 (x 2 ) + y^ j Li j (x)U s+t+ i-j(x), < s e Z. 



With t = 0, this becomes 



s 



T Sj0 (x) = T 0jS (x) - s Li s+ i(x 2 ) + ^2 Lij{x)Li s+1 -.j(x), < s E Z. 
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But for any integers s and t, there holds 



T s Ax) = — — -. r = — 77 r = —T a Ax). 

' n t m s (m — n) ^-^ m s n t [n — m) 

m,n=l v ' m,n=l v ' 



Therefore, 

1 s 

T s,o(s) = - Yl U j( x ) L ^+i-j( x ) - | Li s+1 (x 2 ), < s G Z (8.5) 



On the other hand, 



x <• x 00 x 2n 00 x m-n 00 % n n ~ 1 % r 



T s0 (x) = y-V--— = y — y Y-Y 

n=l m=l n=l m=n+l n=l m=l 



Li^Li^-^E 



rr 14 — ' 7 

n=l 3=1 J 



Comparing this with ()8.5|) gives 



xj n n—1 n—j j s ~ ^ 

J^ — J^— = -Li s+1 (x 2 ) - [Li,(ar) - Li s (x 2 )] Lii(z) - - ^2hi j (x)Li s+1 _ j (x), 

n=l 3=1 J 3=2 

(8.6) 

where in ()8.6|) and what follows, we now require 2 < s G Z because the terms j = 1 and 
j = s in the sum ()8.5|) were separated, and assumed to be distinct. 

Next, note that if n is a positive integer and < x < 1, then 

71-1 

1 - x n = (1 - x) Y. xj < (1 _ x ) n - 

3=0 

Thus, if 2 < s G Z and < x < 1, then 

C^O 71/1 71 \ 00 71 

< [Li s (x) - Li s (x 2 )} Lix(x) = Lix(x) ]T ^ J ^ < (1 - x)Lii(x) ^ 

n=l ra=l 

< (1 - x)\og 2 (l - x). 

Since the latter expression tends to zero in the limit as x — > 1 — , taking the limit in ()8.6|) 
gives 

1 1 s ~ 2 

CO, 1) = - s((s+ 1) - - ]TC(J + 1)C(« - i), 2 < s G Z, 

3=1 
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which is □ 

9. Fourier Series 

The Fourier expansions 

Esin(nt) vr — t , *^ cos(ni) , . . . , . 

-±-± = — and £— ^ = -log|2sin(t/2)| 

n=l n=l 

are both valid in the open interval < t < 2n. Multiplying these together, simplifying, 
and doing a partial fraction decomposition gives 

Esin(ni) 1 1 sin(ni) ^-4 f I 1 \ 1 sin(nt) 
n ^ k 2^ n ^\k n-k) 2 2—* kin - k) 

n=l k=l n=l k=l v 7 n>fc>0 v ' 

1 sin(m + n)t sin(mt) cos(nt) 

2 ^— ' mn ' mn 

m,n=l m,n=l 

"*log|2sin(*/2)|, (9.1) 



(9.2) 



2 

again for < t < 2tt. Integrating (|9.1|) term by term yields 

E ^F 1 E 1 = C(2, 1) + \ / (* - f) log |2 sin(t/2) | dt, 

n=l fe=l ^° 

valid for < 6 < 2vr. Likewise for < 9 < 2tt, 

g oos(n#) = c(3) + /• (0 _ t)log | 2sin( i /2) | ^ (9 3) 

n=l ^ 

Setting 6* = 7r in ()9.2j) and (j9.r>j) produces 

C(2, 1) - C(2, 1) = -\ J\n - t) log |2 sin(t/2)| dt = C(3) ~ C(3) . 
In light of (jl.2j) . this implies 

( (2 1) «3) + C(5) if 1 + (-D- f. 1 1 &) 

n=l m=l v ' 

which is fTKj) . □ 
Applying Parseval's equation to ()9.1|) gives (via [TQ1 112 EZ] ) the integral evaluation 

"27T _00_ -y-y 



1 /^tt ^ rr2 

4- / (--^log 2 (2sin(t/2))dt = ^ I -^ ip = ^C(4). 
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A reason for valuing such integral representations is that they are frequently easier to use 
numerically 



10. Further Generating Functions 

10.1. Hypergeometric Functions. Note that in the notation of (jl.lj) . C(2, 1) is the 
coefficient of xy 2 in 

oo oo oo oo k— 1 / \ 

G(x, y) :=EE x m+1 y n+1 C(rn + 2, {1}") = „ £ x^ £ — J] 1 + V - . (10.1) 

m=0 n=0 m=0 fc=l j=l ^ J ' 

Now recall the notation (y) k := y(y + 1) • • • (y + k — 1) for the rising factorial with 
factors. Thus, 

Substituting this into (jlO.ljl . interchanging order of summation, and summing the result- 
ing geometric series yields the hypergeometric series 

But, Gauss's summation theorem for the hypergeometric function jl, p. 557] |U p. 2] and 
the power series expansion for the logarithmic derivative of the gamma function QjJ p. 
259] imply that 



2^1 



y, -x 
1 — x 



r(i-x)r(i- y ) 



exp{f> fc + ^-(z + #)^H. 

^ k=2 ' 



r(l -x-y) 

Thus, we have derived the generating function equality (see [21] for a g-analog) 
£ $> m+ V +1 C(™ + 2 > «") = 1 - exp (a; fc + / - (x + iU L (i .2) 

m=0 n=0 ^ k=2 ' 

Extracting coefficients of xy 2 from both sides of ()10.2|) yields (jl.2|) . □ 



m-l„,2 



The generalization (jl.HJ) can be similarly derived: extract the coefficient of x 
from both sides of ()10.2j) . In fact, it is easy to see that ()10.2j) provides a formula for 
((m + 2, {l} n ) for all nonnegative integers m and n in terms of sums of products of values 
of the Riemann zeta function at the positive integers. In particular, Markett's formula [SI] 
(cf. also [T2] ) for ((m, 1,1) for positive integers m > 1 is most easily obtained in this way. 
Noting symmetry between x and y in ()10.2J) gives Drinfeld's duality formula 



(( m +2,{in=((«+2,{in (io.3) 
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for non- negative integers m and n, a special case of the more general duality formula (j5.11jl . 
Note that (jl.2j) is just the case m = n = 0. 

Similarly [3H1 2.1b] equating coefficients of xy 2 in Kummer's summation theorem [SHI 
p. 53] |H p. 9] 

t \ = T(l + x/2)T(l+x-y) 
T(l + x)r(l + x/2-y) 



t? i x,y 
2tl[ l + x-y 



yields ffOJl . 



10.2. A Generating Function for Sums. The identity (jl.2|) can also be recovered by 
setting x = in the following result: 

Theorem 7 (Theorem 1 of [TT]). // x any complex number not equal to a positive 
integer, then 

oo 1 n— 1 1 oo 1 

y - y ^L_ = y i . 

2 -— ' n(n — x) 2 -— ' m — x 2 -— ' n 2 (n — x) 

n=l v 7 m=l n=l v y 

Proof. Fix x G C \ Z + . Let 5 denote the left hand side. By partial fractions, 

oo n— 1 / 1 1 



^— ' ^— ' \n(n-m)(m-x) n(n — m)(n — x) 

n=l m=l xv /v ' v /v ' 



n-1 



y 1 y - y - y 1 

^— ' m — x ^— ' n(n — m) ^— ' rim — x) ^— ' n — m 

m=l n=m+l v ' n=l v ' m=l 

oo oo / _. xoo n— 1 

y t y f 1 -i)-y - y~. 

m=l v ' n=m+l x 7 n=l v ' m=l 



Now for fixed m G Z + , 

^— ' \n — m n J n->oo ^— ' \n — m n J ^— ' n iv^oo ^— f TV — n + 1 

n=m+l x 7 n=m+l x 7 n=l n=l 

= ^n' 

n=l 

since m is fixed. Therefore, we have 

oo 1 m 1 oo 1 n— 1 oo * n 1 n— 1 1 

5 = y y y > — - y ( y _ _ y _ 

' m(m — x) ' n ^— ' n(n — x) ' m ^— ' n(n — x) V ' m ^— - ' m 

m=l n=l n=l m=l n=l x m=l m=l 



y , 

2 -— ' n 2 (n — x) 

n=l v 7 



32 



J. M. BORWEIN AND D. M. BRADLEY 



□ 

Theorem [7| is in fact equivalent to the sum formula [SH E2] 

C(«i + 2, a 2 + 1, . . ■ , a r + 1) = C(r + s + 1), (10.4) 

E a i= s 

a;>0 

valid for all integers s > 0, r > 1, and which generalizes Theorem El (|4.1j) to arbitrary 
depth. The identity (jl.2j) is simply the case r = 2, s = 0. A (/-analog of the sum 
formula (jl().4|) is derived as a special case of more general results in j21j. See also [2E]. 

10.3. An Alternating Generating Function. An alternating counterpart to Theo- 
rem [7| is given below. 

Theorem 8. (Theorem 3 of JI]). For all non-integer x 

V (~ 1 ) n f Hn i y x2 \ (-1)" x ^, x _ Ecot^) _ J_\ 

' n 2 — x 2 I n ^— ' n(n 2 — x 2 ) J ' n 2 — x 2 1 2 2x J 

n=l n=l v ' ' n=l v y 



-I n 



o>o odd 



"l^-.c) ^(n 2 -x 2 ) z 



V - x 2 V 1 



,, ,,,„,,. (x* — e*) ' — ', , o (x 2 — o 2 ) 2 
e>o even \ ) o>0 ^ \ > 



Setting x = reproduces f)l .5|) in the form £(2, 1) = Yl^r>o^ n ) 3 - We recor d that 

(-1)" J 7T 

^— ' n 2 — x 2 2x 2 2x sin(7rx) ' 

n=l v 7 

while 

E £S{^> - **) - * «*o»o - = E ^{"» + E^^) } 

n=l k J n=l k n=l v y 

~ ^ (2n- l)((2n- l) 2 - x 2 ) + ^ (n 2 -x 2 ^ 2 ' 

n=l v ' v v y y n=l v 

10.4. The Digamma Function. Define an auxiliary function A by 
xA(x) : = ~x)-\ - x) + 7 ) 2 - |C(2). 

We note, but do not use, that 

1 Z" 00 t (e~ l + e~* (1-a:) ) 1 / r°° e~ l - p-^ 1 " 3 ) \ 2 

* A «4f ( i- e - -r^ dt ) -« 2 '- 
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It is easy to verify that 



tp(l -x) + 7 = ^ 



oo 

X 



n- 

oo 



-f n(x — n) ' 

^ ( ) ^ ( C x _ ^2 ) Z-e 77,2(77, — x) 2 ' ^ 

n=l x 7 n=l v y 



and 



n-l 



^— ' z — ' n(n — a;) ^— ' m 

n=l v ' m=l 



n=0 



Hence, 



n-l 



Now, 



A (a;) = N x N \ 

' n 2 (n — x) ' ri(n — x) ' m(m — x) 

n=l n=l v ' m=l v ' 



00 00 1 n— 1 00 1 n—1 _ 

y i x y I y I = y I Y- 

^— ' r7, 2 (n — x) n(n — x) ^— ' m(m — x) ^— ' n(n — x) ^— ' m 

n=l v ' n=l v ' m=l v ' n=l v ' m=l 



is directly equivalent to Theorem [7| of All). 21 — see [TT] Section 3] — and we have proven 

00 

A(x) = ^C(^ + 2,l)x™, 



n=0 



so that comparing coefficients yields yet another proof of Euler's reduction (|1.3|) . In 
particular, setting x = again produces (jl.2|) . □ 



10.5. The Beta Function. Recall that the beta function is defined for positive real x 
and y by 



B(x,y):= I I" ; il /)•" <H VU)Vl!j) 



r(x + y)' 

We begin with the following easily obtained generating function: 

V- ,nrr M 1 ~ *) 

n=l 
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For m > 2 the Laplace integral (|5.4J) now gives 

(m- 1)! J 1-t 

1 (m-l)log m - 2 (t)log 2 (l-t)- 



2(m - 1)! 7 v ' & v ' & v ' t 
(_1)m &i m - 2) (0), (10.6) 



where 

(cf. ginHD - Since 
d 2 



2(m - 2) 

dy 

B(x, y) = B{x, y) [(if>(y) - ijj{x + y)) 2 + ij/(y) - f(x + y)] 



d 2 



= 2A(-x) 

y=l 



dy 2 
we derive 

6 = m)-^(x+i)) 2 +^'(i)-^(x + i) 

X 

Now observe that from ([10.6)1 . 

C(2,i) = - h(o) = lim (-T-^ + D) 2 _ to ^M^W = - V (1) 

sv ' ; 2 v ; *io 2x xio 2x 2 v v ; 

= C(3). 

Continuing, from the following two identities, cognate to (|10.5|) . 

(-7 - v>(* + 1)) 2 = ( £)(-i) m C(m + 1) * m ) 

^ m=l ' 

oo m— 1 

m=l fc=l 

oo 

C(2) - ij/(x + 1) = E(- 1 ) W+1 (^ + l)C(m + 2) x m , 



we get 



2^(-irCKl)x™- 2 = ^|-||x-- 2 = 6l (x) 



m=2 m=2 

m— 1 



^(-l)" 1 - 1 (m + l)C(m + 2) - ^ C(* + l)C(m " * + 1) 

m=l \ fe=l 



j,m-l 



□ 
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from which Euler's reduction (J1.3J) follows — indeed this is close to Euler's original path. 

Observe that (jl(J.6j) is especially suited to symbolic computation. We also note the 
pleasing identity 

#(x) = ^M (10.7) 

In some informal sense (jl0.7|) generates (jl.3j) , but we have been unable to make this sense 
precise. 

11. A Decomposition Formula of Euler 

For positive integers s and t and distinct non-zero real numbers a and x, the partial 
fraction expansion 

1 v~-v ft + r — l\ 1 r-^ fs + r — 1\ (— l) r 

x s (x-aY~^ ^ ^ \ t-1 ) x s ~ r a t+r + ^ \ s-1 ) a s+r (x - a)^ r 

implies p. 48] [SI] Euler's decomposition formula 

\((t-r)((s + r) 



c( S ,t) = (-i)*x:( t tl7 1 )^--r+r)+x:(-ir( s tl7 1 ) ( 

- (-1)' (' + * ~ 2 ) {C(* + 1) + C(s + 1 - 1, 1)}. (11.2) 

The depth-2 sum formula (|4.1|) is obtained by setting t = 1 in (|11.2|) . If we also set s = 2, 
the identity (|1.2|) results. To derive (|11.2|) from (jll.lj) we follow [ST] , separating the last 
term of each sum on the right hand side of (jll.lj) . obtaining 



x s (x — a) 1 V t — 1 / x s ~ r a t+r ^— ' V s — 1 / a s+r (a; — a0*~ r 



-(-1 



s + t-2\ 1 f 1 1 



s — 1 J a s+t 1 \x — a x y 
Now sum over all integers < a < £ < oo. □ 

Nielsen states (lll.ljl without proof [EH P- 48, eq. (9)]. Markett proves (jll.lj) by induc- 
tion [nH Lemma 3.1], which is the proof technique suggested for the a = 1 case of (jll.lj) 
in [T2J Lemma 1]. However, it is easy to prove (jll.lj) directly by expanding the left hand 
side into partial fractions with the aid of the residue calculus. Alternatively, as in [2*7j 
note that (Jll.lj) is an immediate consequence of applying the partial derivative operator 

oy- 1 i ( d N, - J 



s-l)!\ dx ) (t-l)\\ dy 
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to the identity 

1 1 1 

+ 



xy (x + y)x (x + y)y' 

and then setting y = a — x. This latter observation is extended in [2Z| to establish a 
g-analog of another of Euler's decomposition formulas for ((s,t). 

12. Equating Shuffles and Stuffles 

We begin with an informal argument. By the stuffle multiplication rule P^ll9|l2*^ l2^j 

C(2)C(1) = C(2,1) + C(1,2) + C(3). (12.1) 

On the other hand, the shuffle multiplication rule [151 ITo^ IT§1 |20J |^] gives ab tub = 
labb + bab, whence 

C(2)C(1) = 2C(2,1) + C(1,2). (12.2) 
The identity (jl.2j) now follows immediately on subtracting (j!2.1|) from ()12.2|) . □ 

Of course, this argument needs justification, because it involves cancelling divergent 
series. To make the argument rigorous, we introduce the multiple polylogarithm [UJ H3 
HH]. For real < x < 1 and positive integers s\, . . . ,Sk with x = si = 1 excluded for 
convergence, define 

ni>--->n fc >0 j=l J j=l v r=l l r / A Lsj 

where the integral is over the simplex 

x > 4 1} > • • • > <w > • ■ • > 4 fe) > ■ ■ • > t% > o, 

and is abbreviated by 

"fla'^b, a = f, b=^- f . (12.4) 
i=i 
Then 

C(2)Cx(l) = ^^^y = J] ^+ ^ ^2+^£3' 

n>0 fc>0 n>fc>0 fc>n>0 fc>0 

and 

PX PX PX 

( x (2)^(1) = / aft / 6=/ (2a&6 + &a&) = 2C B (2,l)+C B (l,2). 
Jo Jo Jo 



Subtracting the two equations gives 

[C(2) - C.(2)]C(1) = C-(3) - C(2, 1) + E 



n 2 k 

n>k>0 



GOLDBACH VARIATIONS 



37 



We now take the limit as x — >• 1 — . Uniform convergence implies the right hand side 
tends to £(3) — £(2, 1). That the left hand side tends to zero follows immediately from 
the inequalities 



o<4c(2)-a2)]ai) = ^ / mi-^mi-x 

J X 



dt 
T 



< / log 2 (l -t)dt 

J X 

= (l-a;){l + (l-log(l-a;)) 2 } 



□ 



The alternating case f)l .5|) is actually easier using this approach, since the role of the 
divergent sum ((1) is taken over by the conditionally convergent sum ((1) = — log 2. By 
the stuffle multiplication rule, 



C(2)C(1) = C(2,1) + C(1,2) + C(3), (12.5) 
C(2)C(T) = C(2,T) + C(T,2) + C(3). (12.6) 

On the other hand, the shuffle multiplication rule gives ac lu c = 2ac 2 + cac and ab lu c = 

abc + acb + cab, whence 

C(2)C(T) = 2C(2,1) + C(I,2), (12.7) 
C(2)C(T) = C(2,T) + C(2,T) + C(T,2). (12.8) 

Comparing (|12.5|) with (J12.7)) and (|12.6j) with (|12.8j) yields the two equations 

C(2, 1) = C(T, 2) + 2C(2, 1) - C(T, 2) - C(3), 

C(2,T) = C(T,2)-C(T,2) + C(3). 

Subtracting the latter two equations yields 2£(2, 1) = ((3) + C(3), i-e. ()2.3|) . which was 
shown to be trivially equivalent to (jl.5|) in ^21 D 

13. Conclusion 

There are doubtless other roads to Rome, and as indicated in the introduction we should 
like to learn of them. We finish with the three open questions we are most desirous of 
answers to. 

• A truly combinatorial proof, perhaps of the form considered in |16j . 

• A direct proof that the appropriate line integrals in sections 15.81 and 15.91 evaluate 
to the appropriate multiple of C(3)- 

• A proof of or at least some additional cases of it. 
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